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Abstract 

We calculate the average polar isability of two dimensional mesoscopic rings 
in the presence of an Aharonov-Bohm flux. The screening is taken into ac- 
count self-consistently within a mean- field approximation. We investigate the 
effects of statistical ensemble, finite frequency and disorder. We emphasize ge- 
ometrical effects which make the observation of field dependent polar isability 
much more favourable on rings than on disks or spheres of comparable radius. 
The ratio of the flux dependent to the flux independent part is estimated for 
typical GaAs rings. 

Transport properties of mesoscopic metallic conductors are known for a long time to be 
quite sensitive to the quantum phase coherence of electronic wave functions at low tem- 
perature. More recently it has been shown both theoretically and experimentally [0], that 
the orbital magnetism also presents remarkable signatures of this phase coherence. The 
purpose of the present paper is to discuss to what extend the electrical poiarisability i.e. 
the response of a metallic sample to an electrostatic field is also sensitive to mesoscopic 
effects. This quantity can be experimentally measured by inserting the mesoscopic samples 
into a capacitor. In the linear response regime, the change in the capacitance is directly 
proportional to the average poiarisability a of the particles. If the typical size L of the 
particle is much larger than the Thomas Fermi screening length A s , the poiarisability is 
mostly determined by the geometrical shape of the particle with a small correction of order 
^f- || . If furthermore L < the phase coherence length, one can expect to find a quantum 



correction 5o>q to the polarisability: a = oltf + o~cvq. The typical value of this magnetic field 
dependent quantum correction was estimated from diagrammatic theory || in disks and 
spheres. We have calculated the average of the magneto-polarisability of two dimensional 
mesoscopic rings in the presence of an Aharonov-Bohm flux, with in-plane electric field. In 
this geometry c^tf is of the order of R 3 and the quantum correction is of the order of \ S /W, 
where R and W are the radius and width of the ring. This correction is much larger than 
for disks or spheres of comparable radius with approximatively the same value for a = oltf- 
We consider a time dependent electric field which leads to two different limits whether the 
frequency is smaller or larger than the inverse inelastic time 7. Nevertheless, both u and 
7 are supposed to be small compared to the mean level spacing A. We will consider the 
average value of the polarisability in grand canonical ensemble (GCE) as well as in the 
canonical ensemble (CE). In this letter we present mainly analytical results which have been 
confirmed by numerical simulations, to be published separately [|J . The agreement between 
our findings and the recent results of Efetov on disks and spheres H is only qualitative, (in 
particular, we strongly disagree on the disorder dependence of the effect.) 

We have treated the screening of the applied electric field by the electrons within a mean 
field approximation which leads to the one electron Hamiltonian in an Aharonov-Bohm ring 
(A is the static vector potential): 
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Where $ scr (T > ) = EQF(~r*) is the effective screened potential experienced by the electrons 
in the time dependent external electric field E exp(iu;t) which is assumed to be small enough 
in order that linear response is valid . Vdis is a disordered impurity potential and Vint 
is the mean field potential due to electron-electron interactions. In an Aharonov-Bohm 
geometry $ SC r(~^) and Vi n t(y?) depend on the magnetic flux, as pointed out by Bttiker 0. 
In the present work we neglect Knt("^) an d consider consider single electron eigenstates 
and eigenvalues, but we fully take into account of the screening. <& scr (~r*) is related self 
consistently to the field induced shift of the electronic density 5n(~r*,Eo), according to: 
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The ac electric field induces a dipolar moment which is given by: SP x (t) = 
eTr [(p(t) — Po)X] = aE. Where Tr denotes the trace over all states. X is the projec- 
tion of the position operator along the direction of the electric field and p is the density 
matrix, solution of the master equation in the presence of the time dependent perturbation: 

Op 

i-QT = [H, p] ~ H(P ~ Pequi) (3) 

7 is the inverse relaxation time towards equilibrium. p eq is the instantaneous equilibrium 
density matrix which satisfies: [H, p eq ] = 0. Following the same steps than for the compu- 
tation of the orbital susceptibility of Aharonov Bohm rings, within linear response theory 
the polarisability can be expressed as a function of eigenvalues e a of the unperturbed 
Hamiltonian, matrix elements of the operators X and F{f) taken between the corresponding 
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eigenstates, and populations f a = exp 
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In the case where the dimensions of the sample are at least a few times larger the 

screening length, \ s ~ it is reasonable to describe the screening by a Thomas- Fermi 

approximation, where the induced electronic density is proportional to the potential. For 

a 2D system: XtfC^j r ) = (2?reo/ ^A s )5("r > — r ). The potential in the ring has then the 

form F TF (~r f ) = f(r, A s , W) cos#, where the function f(r) depends on the aspect ratio 

of the ring. For W « R the charge density is nearly equally distributed between the 2 

boundaries of the ring, on a length of the order A s , otherwise the charge is concentrated 

on the external boundary. The Thomas- Fermi polarisability is related to f(r) by: qtf = 
2e 7r 2 
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■ J r 2 f(r)dr. For small A s , we recover the polarisability obtained by solving the Poisson 



equation, whose value depending on the ratio W/R extrapolates between the value of the 
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ID ring am = — — : — - and the value of the disk M a 2 D = eo — R 3 - As, the ratio X s /W 
ln(R/W) 3 

increases, size dependent corrections decrease the value of axF ■ 

In order to estimate the quantum correction to the polarisability let us write the local 
electric susceptibility tensor as following: x{ r u r 2) = Xtf{j'i^'2) + <5x( r i, r 2, </>)• Where 5\ 
is the flux dependent quantum correction. The same for the average screened potential: 
F(~r*,<p) = F TF (Y*) + 8F(l^,(j)). In first order of 8\ and 5F which are self consistently 
related through eq.|2| the quantum correction to the polarisability in the limit (a;, 7 << A) 
can then be expressed as: 

8a(<P) = 2e>5 ( £ \F™( r>)| 2 + £ ^{F^A (5) 



Note that this expression is very similar to the unscreened polarisability |10Hl2| where 
the matrix elements of the operator X are replaced by those of the average, flux independent, 
screened potential. For the discussion of this quantity in the zero frequency limit it will be 
important to distinguish between CE and GCE averages as well as between the order of 
the limits 7 >> to —>■ (thermodynamic polarisability ax) and u » 7 —>■ (dynamical 
polarisability ap). On this point we disagree with the statement of Efetov who claims that 
all the limits should be identical 0. In order to obtain the following results valid at zero 
temperature, we have used the expressions of the CE and GCE average population factors 
< fa~ f/3 > an d < df a /de a > which has been extensively studied in |]8| JT3|JT^1 in the context 
of the ac magnetic response of Aharonov-Bohm rings. We call 5a the flux dependent part 
of a. 

In the GCE: 

' ^ G T° E = S (iTr(F^f) = 
^ = %5 (E^ |F-| 2 ) = -%& (£ a WZ\ 2 ) 
Where A is the average level spacing and Ep is the Fermi energy. In the CE at zero tem- 
perature, the quantities df a /de a = 0, there is no contribution of diagonal matrix elements 
to the polarisability. The thermodynamic and dynamical polarisability are identical. 



5<4* = 5a™ = * 5 I £ \F FF \\l + _*-)) (7) 

Expression (|^) shows that the flux dependent polarisability in the GCE is directly related 
to the matrix elements of the screened potential F TF which can be estimated from the Fourier 
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spatial decomposition: F TF (f) = ^2 F n exp(iq n f) where q n = —u r + —uq and M is 

n=—M/2 

the number of transverse channels in the ring. Using the general relation established by 



McMillan [15] for matrix elements of exp(iqr): 

I f~A 12 A D( l 2 t*\ 

|«p(,^| =- D2q4 + {ea _^ )2 (8) 

D is the diffusion coefficient, for low energy difference the main contribution is given by the 
smallest q vector q for which F = (J f(r)dr)/W = 8R\ s /3tt 2 W ~ 2R/M (taking the 2D 
limit for f(r)). As a result: 

<W™ < 9 > 

E c = hjTD = hD/(2irR) 2 is the Thouless energy and tjj is the diffusion time around the 
ring. This relation allows the generalization of the properties of the X operator established 
iT)| . |T7| . More interesting for our purpose, these matrix elements are also flux dependent. 
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This flux dependence is 0o/2 periodic and presents a minimum at 0o/4, in contrast with 
the average square of the diagonal matrix element which is maximum at 0o/4. The time 
reversal symmetry property of the operator F TF (f) implies that its diagonal matrix elements 
are real, even functions of flux and can be developed in successive powers of cos(2-7r0/0o)- 
1S maximum for multiple values of 0o/2 and minimum for 0o/4 [0o/2] . The opposite 
flux dependence for the non diagonal matrix elements directly follows from the flux indepen- 
dence of the trace Yl a ,p \^I,p\ 2 = HrF{f). Note that flux dependences of the diagonal and 
non diagonal matrix elements of the current operator which is an odd function of flux, are 



just opposite as already noted in fTEfl . In order to estimate quantitatively the flux depen- 
dence of \F FF | 2 , it is useful to relate the diagonal matrix element F FF and the sensitivity of 

eigenvalues e a to the screened electric field: F FF = tt^- Since electric field preserves time 

oEft 
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reversal symmetry [18|], the typical derivative of the energy levels: < |— ^-\ 2 > is propor- 



tional to 1/(3, where the parameter (3 characterizes the symmetry class of the Hamiltonian, 
determined by another parameter (flux, for example), as a result: 

\Fl F a \ 2 {ct) G /A){GUEensemble) = - 1 | 2 (0) (GO E ensemble) (10) 

From equations ^U^JlI] on e deduces the magneto-polarisability in the GCE: 

qg CE (0o/4) - a% CE (0) 8 A A s 8 Ar D X s 
5a/a = a = 3^Y C W = S^^W (U) 

This result implies a positive low field magneto-polarisability in agreement with findings 
of Efetov 0]. We also find an increase of the magneto-polarisability with disorder in the 
diffusive regime, in good agreement with numerical data depicted in figure [I] but at odd 
with reference ||, where a value of the order of (Ap-)|* (r e << r D is the elastic scattering 
time) found for Sa/a. A similar 1/D dependence was also found for the typical value of the 
magneto-polarisability ||. It is also clear from our results that the ring geometry is very 
favourable for the observation of quantum effects on the polarisability. The extrapolation for 
a disk or a sphere of same radius with comparable values of otf yields to 5a/ a = (^g^)^' 
The relative correction is thus reduced by a factor (W/R) 2 for the 2D disk and a factor 
XfW 2 /R 3 for the sphere. We can estimate the relative magneto-polarisability of a ring made 
from a GaAs/GaAlAs hetero junction with the following parameters: L = 8/im, X s = 400 
A, l e = 2/i m , M = 10, E c = 7 A, to be of the order of 310~ 3 . Recent experiments measuring 
the ac complex conductance of such Aharonov-Bohm rings coupled to an electromagnetic 
resonator are sensitive enough to detect their magneto-polarisability 

If we consider now the canonical ensemble (CE), according to eq.|7], the flux dependent 
correction to the polarisability in the CE at T = depends both on matrix elements dis- 
cussed above and energy levels. Assuming that energy levels and matrix elements are indeed 
independent random variables, it is possible to rewrite eq.[7] 
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a GCE + (F TF )2 >a f ^Mfc 

A J e 



(12) 
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where < (F TF ) 2 >a is the average of |-F^| 2 for \e a — £p\ < A and R(e) is the probability 
of finding two levels separated by e whose expression for the different symmetry classes are 
known from random matrix theory [pH[] . It is straight forward to show that: 

j { mp_mp )de= _^_ (13) 

The flux dependence of the average matrix element 5 < (F TF ) 2 >a is of the order of 
——(5 < (Fa, a TF ) 2 > as a result: 

6a CE ~ 5a GCE ^r (14) 

In the limit where g = E c /A » 1, the level repulsion contribution to the magneto- 
polarisability compensates the flux dependence of a ( f ) CE in eqJT^. As a result the flux 
dependence of the polarisability in the CE is strongly reduced compared to the GCE at zero 
temperature. However it is not obvious whether this result obtained with single electron 
wave functions is going to survive if electron-electron interactions are taken into account. 
These differences between CE and GCE averages are also expected to disappear at T > E c . 

We have estimated quantum corrections to the average electrical dynamical polarisability 
of 2D Aharonov-Bohm rings which give rise to a low field positive magneto-polarisability. 
This effect has been shown to be directly related to the flux dependence of the matrix 
elements of the screened potential in the rings. In the CEat low temperature the magneto- 
polarisability is sensitive to level spacing statistics whose contribution cancels the flux de- 
pendence due to matrix elements, it is strongly reduced compared to the GCE value. It is 
possible to evaluate the magneto-polarisability of a typical GaAs/GaAlAs ring to be of the 

&CX 

order of: — = 310~ 3 . This result is encouraging for experiments at least when the limit 

a 

to > 7 is achieved. 
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FIGURES 

FIG. 1. Flux dependent dynamical and thermodynamical ar polarisability in the GCE 
ensemble. These data are obtained from numerical simulations on the Anderson model for a ring 
80 x 8 and 2 different values of disorder: = 1 and = 2 in the diffusive regime. The vertical 
axis has been arbitrarly shifted for the different curves. 
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